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We search for charmless decays of charged B mesons to the three-body final state K^K^n^ . 
Using a data sample of 423.7 fb^^ collected at the T(45) resonance with the BABAR detector, 
corresponding to (465.1 ± 5.1) x 10® BB pairs, we find no significant signal and determine a 90% 
confidence level upper limit on the branching fraction of 5.1 x 10~^. 

PACS numbers: 13.25.Hw, 11.30.Er 
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Charmless decays of B mesons to final states with 
even numbers of strange quarks or antiquaries, such as 
5+ K'^K^tt^ [ij, are suppressed in the standard 
model. Such decays proceed mainly via the b d loop 
(penguin) transition. Hadronic b ^ d penguin transi- 
tions have been observed in the decays B° K'^K^ and 
and their effects have also been seen 
through direct CP violation in charmless B decays, such 
as B° 7r+7r_- and B° vr+Tr-TrOjl, 0]. In con- 

trast to B°-B° mixing, which is a 6 — > d process with 
a change of beauty- flavor quantum number of AF — 2, 
little experimental information exists on AF — lb ^ d 
decay amplitudes. There is still potential for new physics 
effects to be uncovered in these decays. 

The decay B^ ^ K'^K'^tt^ has not yet been observed. 
The upper limit on the branching fraction at 90 % con- 
fidence level (CL) is 3.2 x lO"*^ [i]. A model based 
on the factorization approximation, which makes use of 
heavy-quark and chiral symmetries, predicts a nonres- 
onant branching fraction for B^ ^ K^K'^i:^ of order 
10^^ . Decays via intermediate resonant states can also 
lead to the KIK'^'k'^ final state. This motivates an in- 
clusive analysis incorporating both nonresonant and res- 
onant modes. Based on the measured branching fraction 
B[B+ /2(1270)7r+] = (8.2±2.5) x 10"^ [13, [13, [11 , the 
product branching fraction for 5+ /2(1270)7r+ with 
/a (1270) K^K^s should be around 10"^. Similarly, 
B+ /o(980)7r+ and B+ K*+{&92)K° decays could 
contribute to the K'^K^tt'^ channel. The branching frac- 
tion for B+ K*+(892}K° is predicted to be of order 
10-6 or less [H, [B, [H [H, [13, [Hi • 

Another motivation comes from the recent observa- 
tion of _B+ K^K^TT^ by BaBAR, with an inclusive 
branching fraction of B{B+ K+K-Tr+) = [5.0 ± 
0.5(stat.) ± 0.5(syst.)] x 10"^ [ll|. An unexpected peak 
seen near 1.5 GeV/c^ in the K~^K~ invariant-mass spec- 
trum, which we dub the /x(1500), accounts for approx- 
imately half of the total event rate. If the decay of 
the /x(1500) follows isospin symmetry, then equal rates 
would be expected to K+K' and to K"K°_^ If the 
/x(1500) has even spin, then /x(1500) ^ K°K^ decays 
would resuh in 50% K^^K^ and 50% K°K° final states, 
whereas if the /x(1500) has odd spin, then the K^K*^ fi- 
nal state is forbidden by Bose symmetry. Observation of 
the decay /x(1500) ^ in B+ i^,?i^^7r+ could 

therefore provide information on the spin or the quark 
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content of the /x(1500) and could help to elucidate the 
relationship between this state and similar unexplained 
structures seen in i?+ — > K^K^K^ decays [13, [Hj. 
Structures in the K^K'^ mass spectrum have also been 
observed in two-photon and electron-proton colli- 
sions j23l |. 

We report a search for the decay i?+ K^K^t:'^ . 
The analysis is based on data collected at the PEP-II 
asymmetric-energy e+e^ collider at SLAC. The data 
sample consists of an integrated luminosity of 423.7 fb""'^ 
recorded at the T(4S') resonance (on-peak) and 43.9 fb""'^ 
collected 40 MeV below the resonance (off-peak). The 
on-peak data sample contains (465.1 ±5.1) x 10^ BB 
pairs [25j . 

The BaBar detector is described in detail else- 
where [IHl . Charged particles are detected and their mo- 
menta measured with a five-layer silicon vertex tracker 
(SVT) and a 40-layer drift chamber (DCH) located in- 
side a 1.5 T solenoidal magnet. Surrounding the DCH 
is a detector of internally reflected Cherenkov radiation 
(DIRC), designed for charged particle identification. En- 
ergy deposited by electrons and photons is measured 
by a CsI(Tl) crystal electromagnetic calorimeter (EMC). 
Muons and long-lived neutral hadrons are identified in 
the flux return of the solenoid instrumented with resis- 
tive plate chambers and limited streamer tubes. 

We reconstruct a B^ K^K^tt^ candidate by com- 
bining a pair of mesons and a charged pion. A 
TT+TT^ candidate is formed from a pair of oppo- 
sitely charged tracks with an invariant mass that lies 
within 15 MeV/c^ of the nominal mass [ll|, which 
corresponds to five times the mass resolution. We 
require the ratio of measured lifetime and its un- 
certainty to be greater than 20, the cosine of the angle 
between the line connecting the B and decay vertices 
and the momentum vector to be greater than 0.999, 
and the vertex probability to be greater than 10^®. 
Charged pions coming from the B decay are identified 
with the energy loss (d_B/da;) information from the SVT 
and DCH, and the Cherenkov angle and the number of 
photons measured by the DIRC. The efficiency for pion 
selection is approximately 76 % including geometrical ac- 
ceptance, while the probability for misidentification of 
kaons as pions is less than 15%, up to a momentum of 
4 GeV/c. We require pion candidates not to be consis- 
tent with the electron hypothesis, based on information 
from the clE/Ax, the shower shape in the EMC, and the 
ratio of the shower energy and track momentum. 

Continuum e+e" qq {q — u,d,s,c) events are the 
dominant background. To discriminate this type of event 
from signal we use a boosted decision tree (BDT) 
that combines five discriminating variables. The first of 
these is the ratio of L2 to Lq, with Lj = E^pij cos^ip, 
where 9i is the angle, with respect to the B thrust axis, 
of the track or neutral cluster i, and pi is its momentum. 
The sum excludes the daughters of the B candidate and 
all quantities are calculated in the e~^e~ center-of-mass 
(CM) frame. The other four variables are the absolute 
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value of the cosine of the angle between the B direction 
and the beam (z) axis, the magnitude of the cosine of 
the angle between the B thrust axis and the z axis, the 
product of the B candidate's charge and the flavor of the 
recoiling B as reported by a multivariate tagging algo- 
rithm [28], and the proper time difference between the 
decays of the two B mesons divided by its uncertainty. 
The BDT is trained using off-peak data as well as sim- 
ulated signal events that pass the selection criteria. We 
make a requirement on the BDT output (BDTout) such 
that approximately 96 % of the signal is retained and 
60 % of the continuum background is rejected. 

In addition to BDTout, we distinguish signal from 
background events using two kinematic variables: the 
beam-energy-substituted mass mEs = ^^5/4 — and 
IS.E — Eb — \fs/'^, where y/s is the total e+e^ CM en- 
ergy and [Ebt'Pb) is the four-momentum of the B can- 
didate measured in the CM frame. We select signal can- 
didates that satisfy 5.250 < toes < 5.286 GeV/c^ and 
jAE"! < 0.1 GeV. This region includes a sufficiently large 
range of toes below the signal peak to allow properties 
of the continuum distribution to be determined in the 
maximum likelihood fit. 

Another source of background arises from i3+ 
D^{-^ K2:K2:)t^^ decays, where the final state particles 
are identical to the signal. We reduce this background 
by rejecting any event containing a signal candidate with 
a K^K2; invariant mass in the range 1.82 < Mj^oj^a < 
1.90 GeV/c2. 

The efficiency for signal events to pass the selection 
criteria is 28%, determined with a Monte Carlo (MC) 
simulation in which decays are generated uniformly in 
three-body phase space. We find that approximately 9 % 
of the selected 5+ /\gi4r!?7r+ events contain more than 
one candidate, in which case we choose that with the 
highest i?- vertex probability. We have checked that this 
procedure does not bias the fit variables. In about 2 % 
of the signal events, the B candidate is misreconstructed 
because one of its daughter tracks is replaced by a track 
from the rest of the event. Such events are considered to 
be a part of the signal component. 

We study possible residual backgrounds from BB 
events using MC event samples. These backgrounds arise 
from decays with similar kinematic properties to the sig- 
nal or because particles get lost to, or attached from, the 
rest of the event in the process of reconstruction. The 
BB background modes can conveniently be divided into 
two categories, based on their shapes in toes Siiid AE. 
The first category (BBi) contains only B+ K'^K°K+ 
decays, which peak in toes around the B mass and in AE 
near —0.06 GeV. The second category {BB2) contains 
the remaining BB backgrounds and is mainly combina- 
torial. 

We perform an unbinned extended maximum likeli- 
hood fit to the candidate events using three input vari- 
ables: Toes, ^E, and BDTout- For each category j (sig- 
nal, continuum background, BBi, or BB2), we define a 
probability density function Vj (PDF), and evaluate it 



for each event i: 



(1) 



The signal, continuum background, and BB2 background 
exhibit negligible correlations between toes f^nd AE, and 
so the PDF is further factorized: 

V,{ml,s,AE^) = VMes) ■ V,{AE^) . (2) 
The extended likelihood function is 



^-n^ " n 



(3) 



where nj{nk) is the yield for event category j{k). 

The signal toes distribution is parameterized with the 
sum of a Gaussian and a Crystal Ball function [2^ 
while the AE distribution is parameterized with a mod- 
ified Gaussian function with different widths on each 
side, as well as with additional tails that can be dif- 
ferent on each side. We fix the shape parameters to 
the values obtained from the i?+ K'^K'^tt'^ phase- 
space MC sample. The continuum background toes 
shape is described by an empirical threshold ARGUS 
function, x\/ 1 — exp [— •$(1 — x^)] , with x = 2toes/\/s 
and ^ a free parameter [30| . while the continuum AE 
shape is modeled with a linear function. We describe 
the Toes and AE shapes for the BBi sample with a 
two-dimensional histogram determined from MC events, 
which accounts for correlations between these variables. 
One-dimensional histograms are used to describe the toes 
and AE distributions for the BB2 sample. The BDTout 
distributions for all components are described by one- 
dimensional histograms. These are obtained from MC 
events for signal and the BB background categories. The 
continuum background BDTout shape is determined from 
a combination of off-peak data and on-peak data in a 
continuum-dominated sideband of toes, independent of 
the signal region, from which the expected BB back- 
grounds have been subtracted. 

The free parameters of our fit are the yields of the sig- 
nal, continuum, and two BB background categories, to- 
gether with the ^ parameter of the continuum toes shape 
and the slope of the continuum AE shape. 

We test the fitting procedure by applying it to ensem- 
bles of simulated experiments where events are generated 
from the PDF shapes as described above for all four cat- 
egories of events. We repeat the exercise with qq events 
generated from the PDF while signal events are randomly 
extracted from the MC samples. The BB background 
events are either generated from PDF shapes or drawn 
from MC samples. In all cases, these tests confirm that 
our fit performs as expected. No bias is found for the 
value of the signal yield observed in the data. 

The fit to 16 739 candidate events gives a signal yield 
of 15 ± 15 events, where the error is statistical only. The 
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fit returns yields for the continuum, BBi and -6-62 back- 
ground categories of 15 500 ± 140, 89 ± 25 and 1 140 ± 70 
events, respectively. These are somewhat larger than 
the expected values for the first and last categories and 
smaller for the second, a pattern that can be explained 
by the correlations between these yields. 

The results of the fit are shown in Fig. [TJ In these 
plots the continuum background contribution has been 
suppressed by applying a requirement on the ratio of the 
signal likelihood to the sum of the signal and continuum 
likelihoods, calculated without use of the plotted vari- 
able. The value of this requirement for each plot rejects 
about 97 % of the continuum background while retaining 
63 - 71 % of the signal, depending on the variable. 

We determine the inclusive branching fraction for 
5+ if g/CgTr"*" by dividing the observed signal yield by 
the reconstruction efficiency, the number of BB events in 
the data sample, and the square of the daughter branch- 
ing fraction 6(^2 tt+tt-) = 0.6920 ± 0.0005 [U. 
We assume equal decay rates of T(4S') into B^B^ and 
B°B^ pairs. The value obtained is B{B+ K°K'^7r+) = 
(2.5 ±2.4) X 10^^, where the error is statistical only. The 
statistical sig nificance of the s ignal is 1.1 a, which is cal- 
culated as \/— 2 ln(£o / 'Cmax) , where £max denotes the 
likelihood with the nominal signal yield of 15 events and 
£0 denotes the likelihood with the signal yield fixed at 
zero. 

There is a significant dependence of the selection ef- 
ficiency on the kinematics of the K'^K'^tt'^ final state. 
The nominal efficiency is calculated by assuming a phase- 
space distribution of K'^K^tt'^ events. Since we do not 
know the true distribution, a systematic uncertainty of 
24 % is evaluated from the RMS variation of the efficiency 
across the K'^K'^tt'^ Dalitz plot. Smaller systematic un- 
certainties on the fitted yield arise from uncertainties in 
the PDF shapes (4 events), including possible differences 
between data and MC simulations, which are studied us- 
ing a control sample of B'^ D~{-^ K'^tt~)tt~^ events. 
We assign an uncertainty of 2 events to account for fit 
bias. Other uncertainties on the efficiency arise from 
charged particle reconstruction (0.4%), particle identi- 
fication (1.4%), and the selection (1.8%). The un- 



certainty on the number of BB pairs is 1.1 %. The sys- 
tematic uncertainties are added in quadrature to give a 
total of 38 %. Hence the inclusive branching fraction is 
B{B+ K°K°Tr+) = (2.5 ±2.4 ±0.9) x 10"^ where the 
first (second) error is statistical (systematic). 

Since our result is consistent with no signal, we 
determine a 90 % CL upper limit on the branching 
fraction (Sul)- This limit is calculated by integrat- 
ing the likelihood in the physical region such that 
J^^^ C{x) dxj /q^°° Cix) dx = 0.9, where C{x) is the like- 
lihood function for the signal yield x. We have con- 
firmed that the statistical uncertainties from the fit are 
Gaussian, to a good approximation. We therefore as- 
sume a Gaussian behavior for the overall likelihood, with 
a width calculated from the sum in quadrature of the 
statistical and systematic uncertainties. Our result is 
B{B+ K°K^Tr+) < 5.1 x lO^'^ at 90% CL. 

The lack of signal in this decay mode contrasts with 
that observed for B+ K+K-it+ This result dis- 
favors models in which the /x(1500) has even spin and 
decays with isospin symmetry. If the /x(1500) is con- 
firmed to have even spin in future measurements, this 
may indicate a non-gg nature of this state. 

In conclusion, with a data sample of 423.7 fb~^, we 
have performed a search for the decay B^ — > K^K^-k^ . 
We observe no significant signal and set a 90 % confidence 
level upper limit on the branching fraction of 5.1 x 10~^. 
This result provides useful information for the under- 
standing of low energy spectroscopy. 
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